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We investigate Josephson and proximity effects on the surface of a topological insulator on which 
superconductors and a ferromagnet are deposited. The superconducting regions are described by 
the conventional BCS Hamiltonian, rather than the superconducting Dirac Hamiltonian. Junction 
interfaces are assumed to be dirty. We obtain analytical expressions of the Josephson current and 
the proximity-induced anomalous Green's function on the topological insulator. The dependence 
of the Josephson effect on the junction length, the temperature, the chemical potential and the 
magnetization is discussed. It is also shown that the proximity-induced pairing on the surface of 
a topological insulator includes even and odd frequency triplet pairings as well as a conventional 
s-wave one. 
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I. INTRODUCTION 

Topological insulator offers a new state of matter topo- 
logically different from the conventional band insula- 
tor Edge channels or surface states of the topolog- 
ical insulator are topologically protected and described 
by Dirac fermions at low energies. The nature of the 
surface Dirac fermion of the topological insulator mani- 
fests itself in interesting phenomena such as the quantized 
magneto-electric effect^, giant spin rotation 9 -, magnetic 
properties of the surface stated, magnetization dynam- 
ics 1 -^—, magneto-transport phenomena^ 4 - - — , and super- 
conducting proximity effect^Sr— . 

There have been a great and increasing interest on 
topological insulators attached to superconductors. In 
particular, Majorana fermions emerging in these systems 
have been intensively investigated. 25,26 When supercon- 
ductor/ferromagnet junctions are deposited on the topo- 
logical insulators, surface Dirac fermions acquire a do- 
main wall structure of the mass. At the domain wall, Ma- 
jorana fermions emerge as a zero energy bound stat o 27 ' 28 . 
Majorana fermions have received much interest from the 
viewpoint of fundamental physics and also fault-tolerant 
quantum computing due to their exotic propertie a 25 ' 26 . 
It has been also shown that Majorana bound states 
crucially influence the Josephson effect. The current- 
phase relation shows 4tt periodicity, i.e., sin(</>/2) with 
the phase difference across the junction q&-~—. In previ- 
ous works,— -26 ' 28 it is assumed that the Dirac fermions 
become superconducting due to the proximity effect, and 
the junctions between superconducting and normal Dirac 
fermions are considered. In this paper, we take a different 
modeling of the same system. We consider the coupling 
between conventional superconductors and a topological 
insulator, rather than that between superconducting and 
normal (or magnetic) Dirac fermions^Ir— Namely, tun- 
neling between the Schrodinger electrons and the Dirac 
fermions is explicitly taken into account. Here, the su- 
perconductors are topologically trivial and hence, in this 
setup, there appear no Majorana fermions. 

In this paper, we study Josephson and proximity ef- 



fects on the surface of a topological insulator on which 
superconductors and a ferromagnet arc deposited. The 
superconducting regions are described by the conven- 
tional BCS Hamiltonian, rather than superconducting 
Dirac electrons. We consider disordered junction inter- 
faces in contrast to the previous works ^ r 26 i 28 We obtain 
analytical expressions of the Josephson current and the 
proximity-induced anomalous Green's function on the 
topological insulator. The dependence of the Joseph- 
son effect on the junction length, the temperature, the 
chemical potential and the magnetization is discussed. It 
is also shown that the proximity-induced pairing on the 
surface of a topological insulator includes even and odd 
frequency triplet pairing o 34 ' 35 as well as a conventional 
s-wave one. 

Previous works on the Josephson effect on the surface 
of a topological insulator are mostly based on the ap- 
proach which takes into account only the contribution 
from the Andreev bound states! 23 ' 24 ' 28 This holds for 
short junctions with d <C £, where d and £ are the junc- 
tion length and the superconducting coherence length, 
respectively. 36 In this paper, we adopt the functional in- 
tegral method 37 ' 38 which is applicable to any length of 
the junctio n 39 ' 40 , thus allowing to study the asymptotic 
behavior of the Josephson current for d — > oo. 



II. FORMULATION 

We consider superconductor/topological insula- 
tor/superconductor junctions where a ferromagnet is 
also attached to the topological insulator, as shown 
in Fig. [1] (a). Junctions without the ferromagnetic 
region (Fig. Q](b)) can be considered just by setting the 
exchange field to zero in the ferromagnetic region. The 
total Hamiltonian of the system reads 



H — Hr, + Hfj + Hm + H[ 
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FIG. 1: (Color online) Schematic of the model. 
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Here, A and 

<Pl(r) are the magnitude of the gap function and the 
phase of the left (right) superconductor, respectively. 
Also, m is the exchange field, and a and r are Pauli 
matrices in spin and Nambu spaces, respectively. #l(.r) 
represents the Hamiltonian on the left (right) supercon- 
dutor, while Hm is the Dirac Hamiltonian with the ex- 
change field. Note that the superconductors are de- 
scribed by the Schrodinger electrons and topologically 
trivial. Hence, in this setup, no Majorana fermions 
emerge. — Ht is the tunneling Hamiltonian between the 
superconductors and the surface of the topological insu- 
lator which is treated as a perturbation, ^l(r) is the 
position of the interface between the left (right) super- 
conductor and the topological insulator. We consider 
the incoherent tunneling model where the spin is con- 
served but the momentum is not conserved upon tun- 
neling at the interface. This modeling is applicable to 
junctions with imperfect dirty insulating barriers. In real 



space representation, the tunneling matrix element reads 
t5(r — r^(^)). The average of the position vectors is as- 
sumed to give (\rjt — r^]) = d4£ The calculated results 
are averaged over the positions of and at the in- 
terfaces. 

The partition function is then given by 



DipDip exp 



-EH- 

{k} 



(5) 



where ip = (4>k L ,4>k,4>-k R )- Gq is the bulk Green's func- 
tion while T is a tunneling matrix. See the Appendix for 
their explicit forms. The free energy of the system can be 
calculated as F = —T\nZ where T is the temperature 
of the system. The leading contribution to the Joseph- 
son current is given by the fourth order with respect to 
the tunneling Hamiltonian (see the Appendix for the de- 
tails of the calculation) . The Josephson current is then 
calculated as 
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where v, V, u n , and K v {z){y = 0,1) are, respectively, the density of states at the Fermi level, the area of the 
surface of the topological insulator sandwiched between the superconductors, the fcrmionic Matsubara frequency, and 
the modified Bessel function. Also, k F is defined by hv F k F = y (w n — i[i) 2 + m 2 and the branch is taken so that 
Refc^ > 0. Here, cp = cpn — <pi, is the phase difference across the junction. It is seen that the Josephson effect is 
independent of m x , and m y shifts the phase difference ! 23 ' 24 
The critical current Iq can be written as 
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FIG. 2: Critical Josephson current as a function of m z for 
T/T c = 0.1, d/i = 1 and fi/T c = 100. 



where Tc is the superconducting transition temperature 
and 
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III. RESULTS 
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A. Josephson effect 

In what follows, we will study the critical Josephson 
current using Eq.©. We consider a temperature depen- 
dence of the gap of the BCS type modeled by^ 2 - 



A(T) = A(0) tanh ( l.Uy/Tc/T-1 



(9) 



1. The effect of the exchange field 



Here, let us study the effect of the exchange field. 
As seen from Eq.©, the Josephson effect is indepen- 
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FIG. 3: Critical Josephson current as a function of temper- 
ature of the system with m z — and n/Tc = 100 for (a) 
d/i = 0.1 and (b) d/£ = 5. 



dent of m x , and m y shifts the phase difference. Since 
the inplane exchange field corresponds to the shift of the 
momentum^ the effect of the inplane exchange field can 
be reduced to the phase factor (which can be seen by 
proper transformations in Eq. (|A8[l ) and hence we find 
the phase shift proportional to m y . 

To see the effect of the ^-component of the exchange 
field, m z , we plot the dependence of the critical Joseph- 
son current on m z in Fig. [2] for T/T c = 0.1, d/£ = 1 
and n/T c — 100 where £ = Tiv F /T c is the superconduct- 
ing coherence length. With increasing m z , Iq increases 
and for m z > fj, 7 the current is strongly suppressed. This 
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FIG. 4: (Color online) Critical Josephson current as a func- 
tion of the distance between the superconductors d for m z — 
0, n/Tc = 100 and several temperatures. 
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FIG. 5: (Color online) Critical Josephson current as a func- 
tion of the chemical potential /i for m z — 0, T/Tc = 0.1 and 
several d. 



is because for m z > fi, the Fermi level lies inside the 
mass gap and hence the surface state of the topological 
insulator is no more metallic. Note that the the surface 
spectrum is given by Em = ± \Jk 2 + m 2 ± fi (here we 
set m x — m y = for simplicity). Also, it should be 
noted that since the superconductors are topologically 
trivial, there is no edge state, or Majorana bound state. 
This is in stark contrast to similar junctions where su- 
perconducting regions are described by Dirac fermions: 
there, the Majorana bound state can carry the Josephson 
current 



m = 0. This corresponds to the junction illustrated in 

Fig. ni (b). 

Figure |3] depicts T dependence of the critical Joseph- 
son current with /i/Tc = 100 for (a) d/£ = 0.1 and (b) 
d/^ = 5. For short normal segment d/£ — 0.1, the behav- 
ior is similar to that of the conventional Josephson junc- 
tions through an insulating barrier, i.e., tanh(A/2T) 43 . 
For large d, the critical current shows an exponential 
decay. This can be also obtained as follows. Using 
the asymptotic form of the modified Bessel function for 
\z\ » 1: 
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Notice that the n — component in the Matsubara fre- 
quencies has a dominant contribution to the Josephson 
current. This form shows a typical exponential decay 
of the critical Josephson current for Td/hvF ^S> 1: the 
asymptotic behavior of the Josephson current (the expo- 
nential decay) has the same form as that governed by the 
Schrodinger electrons! 40 ' 44 

In Fig. [U we show d dependence of the critical Joseph- 
son current for [a/Tc = 100 and several temperatures. 
For large d or at high temperature, we see an exponen- 
tial decay of the critical Josephson current. This is also 
consistent with the above analytical expression. 

Figure O shows /i dependence of the critical Josephson 
current for T/Tq = 0.1 and several d. It is found that 
with the increase of /i, the critical current monotonically 
decreases. This is because the proximity effect is sup- 
pressed by increasing the chemical potential /i as will be 
shown in Eq. (|T7|) . For large d, the Josephson current is 
inversely proportional to the chemical potential as seen 
from Eq. (jlip. Experimentally, the chemical potential 
can be tuned by chemical doping 46 or gating^. 

Recently, Josephson supercurrent through a topologi- 
cal insulator surface state has been observed^ The de- 
pendence on T and d shown in Figs. [3] and HI is qualita- 
tively consistent with the experimental data. A quanti- 
tative difference would come from the fact that the bulk 
states of the topological insulator also contribute to the 
Josephson current because the chemical potential of the 
sample used in Ref. |45| probably crosses the bulk bands. 



B. Proximity effect 



2. Superconductor /topological insulator /superconductor- 
junction 

Now, we will investigate the Josephson junction char- 
acteristics. Since the effects of the exchange field have 
been clarified, here let us consider the junction with 



In this subsection, we will investigate the proximity 
effect in a topological insulator/s-wave superconductor 
junction. Proximit y eff ect in this junction has been inves- 
tigated in Refs. l3ll - [33l The tunneling between the super- 
conductor and a bulk topological insulator has been con- 
sidered, and the validity of the Fu-Kane model has been 
discussed, based on mostly numerical approaches^ - — 
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Here, we consider the tunneling between the supercon- 
ductor and the surface of the topological insulator. Us- 
ing the functional integral method, we derive analytical 
expressions of the proximity-induced anomalous Green's 
functions on the topological insulator. 

Now, let us consider superconductor/topological insu- 
lator bilayer junctions. To do so, let us remove the degree 
of freedom of the right superconductor from the above 
formulation (see also the Appendix). Without loss of 
generality, we can set (pi — and = 0. The Green 
functions can be calculated as 



where if) 
we have 



Performing the functional integral, 



G = (Go" 1 - T)~ 



(13) 
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where 



G=[ G ' L a> )^=( G t ° ),f= ° Tl2 ). (14) 



G' M J^° V G mJ \T 21 

The leading contribution is given by the second order with respect to the tunneling matrix. The anomalous Green's 
function on the surface of the topological insulator f' M in the second order in t can be represented as 

f'u = - t2 9M(Ku n )g M (k,uj n )^2f L (k L ,uj n ) (15) 

kt 

TTvAt 2 
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(16) 



with kj_ = (k y , —k x ,0). Here, spin-singlet pairing is proportional to the unit matrix in spin space while spin-triplet 
pairing is proportional to the Pauli matrix er. Therefore, it is seen that both singlet and triplet pairings are induced 
on the surface of the topological insulator. The generation of the triplet pairing reflects the symmetry breaking in 
spin space. 48 In particular, for m = 0, we have 



_ irv At 2 - [ujj + pi 2 + {Tiv F k) 2 } - 2/^Fk_L ■ cr 
hl ~ + A 2 K + M 2 ) 2 + {hv F kf + 2(Hv F k) 2 (u; 2 n - f i 2 ) ■ ' ' ' 



We see that in the limit of \i — > oo, we have f' M — > 0. This explains the suppression of the Josephson curernt with /i 
in Fig. [SJ It is also found that even in the absence of the exchange field, triplet pairing is induced on the surface if 
/i ^ 0, which is consistent with Refs. I3lll32l (see also Ref. [49t ). In previous works, it is assumed that by attaching an 
s-wave superconductor to a topological insulator, the same s-wave superconductivity is induced on the surface 2 ^ -26 ' 28 . 
Here, we find that not only s-wave singlet superconductivity but, in general, triplet p-wave superconductivity is also 
induced on the surface of the topological insulator . 31 ' 32 Also, we here assume clean surface states on the topological 
insulator. If the surface is in the diffusive regime, it is expected that odd frequency triplet s-wave superconductivity 
is induced on the topological insulator^ 

Let us focus on the case with fi — but finite exchange field. The anomalous Green's function then becomes 

f' M = F{k x , k y ,uj n ) [-cj 2 - (hv F k) 2 + m 2 + 2icu n m ■ a + 2ihv F (kx x m) • cr] , (18) 

., , . TlvAt 2 

F(k x ,ky,uj n ) 



[w 2 + [hv F ) 2 {{k y + m x ) 2 + (k x - m y ) 2 } + m 2 } [w 2 + {Tiv F ) 2 {{k y - m x ) 2 + (k x + m y ) 2 } + m 2 



(19) 



Note that F(k x ,k v ,u> n ) is an even function of k(= (^,^,0)) and w n . We find that the component propor- 
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tional to — — (hvpk) 2 + m 2 represents a singlet s-wave 
superconductivity while that proportional to 2iw„m • a 
is triplet and odd in uj n , namely odd frequency triplet 
s-wave pairing2^2S,. The component proportional to 
2i(k± x m) • a corresponds to a triplet p-wave super- 
conductivity. 

IV. SUMMARY 

In this paper, we have investigated Josephson and 
proximity effects on the surface of a topological insula- 
tor on which superconductors and a ferromagnet are de- 
posited. We have described the superconducting regions 
by the conventional BCS Hamiltonian, rather than the 
superconducting Dirac Hamiltonian. We have presented 



analytical expressions of the Josephson current and the 
proximity-induced anomalous Green's function on the 
topological insulator. The dependence of the Josephson 
effect on the junction length, the temperature, the chem- 
ical potential and the magnetization has been discussed. 
It has been also shown that the proximity-induced pair- 
ing on the surface of a topological insulator includes even 
and odd frequency triplet pairings as well as a conven- 
tional s-wave one. 

This work was supported by Grant-in-Aid for Young 
Scientists (B) (No. 23740236) and the "Topological 
Quantum Phenomena" (No. 23103505) Grant-in Aid for 
Scientific Research on Innovative Areas from the Ministry 
of Education, Culture, Sports, Science and Technology 
(MEXT) of Japan. 



I 

Appendix A: Calculation of the free energy 



Here, we present the details of the calculation of the free energy of the junctions. The unperturbed Green's function 
Go is represented by a 12x12 matrix as 

Gl 1 
Go" 1 = ( Gl] 1 ,(A1) 
G^ 1 



G 
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V JL(R) 9l(R) 

= (9 M 
9 M 



By performing the functional integral, we have the free energy of the junctions of the form 



F = -T\nZ = -TTrln 
where T is the tunneling matrix given by 



Tia 
T= | T 2 i T 23 
T 32 



,(A2) 
.(A3) 

(A4) 
(A5) 



with T12 = te l ( k - ki > ri cr ® t 3 = and T 23 = te^ R -^- rR a <g> t 3 = T 3 * 2 . 

The leading contribution is given by the fourth order with respect to the tunneling element, which is calculated as 

(g t) 4 

= -TTr ^2 G L (k L ,w„)Ti 2 G M (k, uj n )T23G R (k R ,Lo n )T 32 G M (k' , u n )T 21 
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Here, we have used the following relations 
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for a > and Kek > 0, where Jq(x) is the Bessel function. 
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